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Abstract 



y-^ I We have studied radiative B meson decays into higher kaonic resonances 

(^ ' decaying into a two-body or three-body final state, using a data sample of 21.3 

L^ ■ fb~ recorded at the T{4S) resonance with the Belle detector at KEKB. For 

the two-body final state, we extract the B -^ i^|(1430)7 component from an 

analysis of the helicity angle distribution, and obtain B{B^ — s- /f|(1430)'^7) = 

^ \ (1.26 ±0.66 ±0.10) X 10"^. For the three-body final state, we observe a B —^ 

I ' Kmr'y signal that is consistent with a mixture of B ^ K*tt^ and B -^ Kpj. 

bj ■ This is the first time that B — > K*Tr^ and B —>■ Kpj have been observed 

r^ • separately. We find their branching fractions to be B{B -^ K*tt^; Mk*-k < 

2.0 GeV/c2) = (5.6±1.1±0.9) x 10"^ and B{B -^ Kpr, Mxp < 2.0 GeV/c^) = 

V ! (6.5 ± 1.7+|-2) X 10-^ respectively. 
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I. INTRODUCTION 

Radiative B meson decay through the b —>■ s'j process has been one of the most sensitive 
probes of new physics beyond the Standard Model (SM). The inclusive picture of the b —>■ s'j 
process is well established; however, our knowledge of the exclusive final states in radiative 
B meson decays is rather limited. To date, we know that around 15% of 6 ^ 57 can 
be accounted for by S ^ i^* (892)7 decays. In addition, a relativistic form- factor model 
calculation |1|] predicts that another 20% of the 6 — > 57 process should hadronize as one 
of the seven known higher kaonic resonances (Table |I|). CLEO has already reported an 
indication of the B -^ i^g (1430)7 signal [Q. Precision measurement of the inclusive b ^ s'j 
branching fraction will require detailed knowledge of such resonances, for example to model 
the decay processes into multi-particle final states. In this analysis, we study radiative B 
meson decay processes into higher kaonic resonances, which subsequently decay into two- 
body or three-body final states. 

We have analyzed a data sample that contains 22.8 x 10® BB events. The data sample 
corresponds to an integrated luminosity of 21.3 fb^^ collected at the T(4S') resonance with 
the Belle detector |^ at the KEKB e~^e~ collider p. The beam energies are 3.5 GeV for 
positrons and 8 GeV for electrons. 

Belle is a general purpose detector with a typical laboratory polar angular coverage be- 
tween 17° to 150°. Charged tracks are reconstructed with a 50 layer central drift chamber 
(CDC), and are then extrapolated and refitted with a three layer double sided silicon vertex 
detector (SVD) to provide precision track information for the decay vertex reconstruction. 
Particle identification, namely discrimination of kaons from pions, is provided by combining 
information from silica aerogel Cherenkov counters (ACC) and a time-of-flight counter sys- 
tem (TOF), together with specific ionization {dE/dx) measurements from the CDC. Photons 
are measured with an electromagnetic calorimeter (ECL) of 8736 CsI(Tl) crystals. These 
detectors are surrounded by a 1.5 T superconducting solenoid coil. 

TABLE I. Predicted branching fractions for radiative B decays into kaonic resonances. 

Theoretical prediction [xlO~^] 
mode Veseli-Olsson |l| Ali-Ohl-Mannel S sub-decay modes 

Kir [~100%] 

Kp [42%], i^o*(1430)7r [28%], /^*(892)7r [16%] 
i^*(892)^ [94%], Kp [3%] 
i^*(892)^ [>40%], Ktt [6.6%] 
Ktt [49.9%], K*{892)tt [24.7%] 
iv:*(892)7r 
Kmr, K<f> 
Ktt [38.7%], Kp [31.4%], K*{892)tt [29.9%] 



B- 


-^ K* (892)7 


4.71 ± 1.79 


1.4- 4.9 


B - 


-^Ki (1270)7 


1.20 ±0.44 


1.8- 4.0 


B - 


-^Ki (1400)7 


0.58 ±0.26 


2.4- 5.2 


B - 


-^ K*{UW)'y 


1.14 ±0.18 


2.9- 4.2 


B - 


-> i^2* (1430)7 


1.73 ±0.80 


6.9 - 14.8 


B - 


-^ K2 (1580)7 


0.46 ±0.11 


1.8- 2.6 


B - 


-^Ki (1650)7 


0.47 ±0.16 


(not given) 


B - 


-^K*{l680)-f 


0.15 ±0.04 


0.4- 0.6 



II. EVENT RECONSTRUCTION 

We select events that contain a high energy (1.8 to 3.4 GeV in the T(45') rest frame) 
photon (7) candidate inside the acceptance of the barrel ECL (33° < 6^ < 128°). The 
photon candidate is required to be consistent with an isolated electromagnetic shower, i.e., 
95% of its energy is concentrated in the central 3x3 crystals and there is no associated 
charged track. We combine it with other photon clusters in the event and reject it if the 
invariant mass of the pair is consistent with a 7r° or 77. 

Kaonic resonance {Kx) candidates are formed by combining one kaon with one or two 
pions, in the two-body K^ti^ , K^tx^ and K^ti^ final states and in the three-body K^-k^-k^ 
final state. (Here and throughout the paper, charge conjugate modes are implicitly included.) 
For every charged particle track of good quality, the particle identification information is 
examined. A likelihood ratio for the kaon and pion probabilities is calculated by combining 
information from the ACC, TOF and dE/dx systems. We apply a tight cut with an efficiency 
of 85% for charged kaon candidates and a loose cut with an efficiency of 97% for charged 
pion candidates. Neutral pion (tt^) candidates are reconstructed from pairs of photons 
that satisfy the following requirements: the invariant two photon mass is consistent with 
a 7r°, each photon has more than 50 MeV energy, the opening angle of two photons is 
less than 40° and one of the photons should deposit its energy in more than one crystal. 
The 7r° momentum is recalculated with a 7r° mass constraint. Neutral kaon (K^ -^ tt+tt^) 
candidates are reconstructed from two oppositely charged tracks, whose invariant mass is 
consistent with K^. We require that the K^ candidate form a vertex displaced from the 
interaction point and lie in a direction consistent with the K^ momentum. 

We reconstruct B meson candidates by forming two independent kinematic variables: 
the beam constrained mass and the energy difference. Both variables are calculated in the 



T(4S') rest frame. The beam constrained mass is defined as Mbc = \/(-Ebeam)^ — \Pkx + P7P) 
in which the photon energy is constrained to be E^ = -Ebeam — ^Kx- This constraint 
improves the Mbc resolution by about 20%, resulting in an Mbc resolution of 3 MeV/c^. The 
energy difference, AE = E^x + E^ — -Ebeam, has an asymmetric resolution, mainly due to 
energy leakage from the counters and energy loss in the inner material. We apply a cut of 
— 100 MeV < AE < 75 MeV, which is about a 2a cut on the higher side, and which removes 
around 25% of events on the lower side. 

The largest background source is continuum light quark-pair {qq) production, in which 
the high energy photon mainly comes from the initial state radiation (e^e^ — *■ qqj) and 
decays of neutral hadrons (vr^, rj, . . .). In order to reduce such background, we form a Fisher 
discriminant which we call the Super Fox- Wolfram (SFW) variable 0, 

4 

SFW = a2RT'°" + a^RT'°" + E A^r"", 

1=1 



^major ^ Y.i\Pi\\P-,\Pl{cOS Oj^) ^^^^^^ ^ Ei,j \Pi\\Pj\PlicOSe, 
^ EibilKI ' ^ E^,j\P^\\Pj\ 



where / runs from 1 to 4 for the Legendre function Pi and i, j run over all the neutral 
and charged tracks that are not used to form the B candidate, and coefficients a^ and (5i 
are optimized to maximize the discrimination. The variable is calculated in the signal B 



candidate rest frame rather than the T(4S') rest frame in order to ehminate any correlation 
with Mbc- For further background suppression, we form a hkehhood ratio (LR) from the 
SFW variable and the B meson flight direction (cos^^), 

Ln{bt W , cosUb) 



psis(SFW)p'^cos 9b) + p^s(5Fiy)pbg(cos Ob) ' 

where p^'^ and p°^ are the probability density functions (PDF) for signal and background. 
The PDFs for the SFW are parametrized with an asymmetric Gaussian function. For cos6b, 
we use p^^^{cos9b) oc sin^6'B for signal and a flat distribution for background. We require 
LR > 0.7 for the two-body Kx final states and LR > 0.9 for the three-body Kx final states. 
The LR cut efficiency is determined from the B -^ Dtt data samples using parametrizations 
determined from the signal and continuum qq Monte Carlo (MC) simulation samples. 

The background from other B meson decays is examined with the corresponding MC 
samples. We find a negligible contribution from hadronic charmless decays. Background 
from 6 — i> c decays makes a non-negligible contribution in the sideband of negative AE espe- 
cially at high Kx mass {M^x > 1-5 GeV/c^), but does not contribute in the signal region. 
For the three-body final states, there is a small contribution from B -^ K*{892)'j decays 
especially in the positive AE region; this contribution is removed by rejecting candidates if 
Mbc and AE calculated from K~^TT~'~f falls into the K*'~f signal region. Cross-feed from other 
b ^ s'-f final states is not negligible especially for the Ktttt'j final state. These contributions 
are estimated using an inclusive b —>■ sj MC sample, and subtracted from the signal yield. 

We extract the signal yield from a fit to the Mbc distribution. The shape is modeled 
as a sum of a Gaussian function for the signal and a threshold-type function (ARGUS 
function 0) for the combinatoric background contribution. The normalizations are fioated 
for both components. The signal function is fixed from the B -^ Dn data. The background 
function is determined from AE sideband data in the range 0.1 GeV < AE < 0.5 GeV. 
Using a continuum qq MC sample, we check that the background shape has no visible 
correlation with AE. In order to estimate the systematic error of the fitting procedure, we 
vary the mean and width of the signal shape by ±lcr and use the background shape from 
other sources, namely the continuum MC sample and the LR sideband {LR < 0.3) in which 
the signal contribution is negligible. We assign the largest deviation as the systematic error 
of the signal yield. As a cross-check, we fit the AE distribution with a signal shape from 
MC and a linear function for background with a slope determined from the Mbc sideband 
(5.2 GeV/c^ < Mbc < 5.26 GeV/c^), and obtain consistent results. 

III. ANALYSIS OF B ^ K| (1430)7 

The B — >■ 7^2(1430)7 analysis is performed by requiring the Kir invariant mass to be 
within ±125 MeV/c^ of the nominal 7^2(1^30) value. The results of fits to the Mbc distri- 
butions are shown in Fig. |I|, separately for the neutral and charged modes. The signal yield 
is 29.1 ± 6.7lil'g events for the neutral mode, of which the contribution from other b ^ S'y 
decays is estimated to be 0.4 ± 0.3 events; we also see some indication of a signal in the 
charged mode. These are consistent with the yields calculated from the AE distributions 
shown in Fig. 0. 
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FIG. 1. The Mbc distributions for B° -^ i^|(1430)°7 (left) and B+ -^ i^|(1430)+7 (right) 
candidates. The soUd Une is the fitting result. The background component is shown as the dashed 
line. 
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The event selection efficiency is determined from a MC sample that is calibrated with 
high statistics control data samples. Table |I| summarizes the sources of systematic error; 
a detailed description is given in Ref. [|]. The signal efficiency is (6.99 ± 0.55)% for B^ -^ 
ii'2(1430)°7, including the sub-decay branching fractions. 

In order to distinguish the B — ^ ^^2 (1430)7 signal from B -^ K* (1410)7 ^^^ non- 
resonant decays, we examine the helicity angle distribution for the signal candidates. All 
three modes have different helicity distributions: cos^ 6'hei — cos^ 6*1161 for K2 (14:30), 1 — cos^ 6'hei 
for i^*(1410) and uniform for non- resonant decay. We divide cos6'hoi into 5 bins, and extract 
the yield from fits to the Mbc distribution for each bin (Fig. ^. This distribution clearly 
favors B — > ii'2 (1430)7. We fit the cos^hci distribution with a function which is artificially 
parametrized to avoid a negative yield of i^*(1410) or non-resonant component, and obtain 
20.1 ± 10.5 events for the B —>■ fC|(1430)7 component. After subtracting other b ^ s^ 
contributions, this leads to a 5'^ ^ i^|(1430)°7 branching fraction of 



B{B^ -^ is:*(1430)°7) = (1.26 ± 0.66 ± 0.10) x 10" 



TABLE II. Contributions to the systematic error. 



B -^ i^2* (1430)7 



B -^ Kmrj 



photon reconstruction 
charged track reconstruction 
charged kaon selection 
charged pion selection 
L.R. + tt^/t] veto + vertex 
sub branching ratio uncertainty 



5.3% 
3.0% 
1.7% 
0.6% 
2.4% 
2.4% 



5.3% 
4.4% 
1.7% 
1.2% 
5.1% 



total 



7.2% 



9.1% 
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FIG. 3. The background subtracted 
i^|(1430) helicity angle distribution. The solid 
curve is the fitting result. The theoretical curve 
for B -^ i^l (1430)7 {B -^ iC* (1410)7) is shown 
as the dashed (dotted) line. 
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This result agrees with the prediction from the relativistic form-factor model of Veseli and 
Olsson ||I|], but is much lower than that from the non- relativistic form- factor model of Ali, 
Ohl and Manuel |]. We also obtain the upper limit B{B° -^ K*{UlOf-f) < 8.03 x 
10^^ (90% CL.) by conservatively neglecting the non-resonant component in the fitting 
procedure. 

The background subtracted Kn invariant mass distribution for B -^ Ktt'j is obtained 
by a similar method. We divide the Mktt spectrum into 100 MeV/c^ bins, and extract the 
signal yield from the Mbc distribution for each bin. In Fig. ^. we see a clear enhancement 
around 1.4 GeV/c^, which supports the conclusion that the B -^ 7^1(1430)7 contribution 
dominates. 



IV. ANALYSIS OF B ^ Kxj -^ K-rttj 

The selection criteria used to reconstruct the B -^ Krni'j decay are identical to those used 
in the analysis oi B ^ i^2 (1430)7, unless explicitly stated otherwise. The Kx candidate 
is reconstructed from K^ti^ti^, and required to have a mass between 1.0 GeV/c^ and 2.0 
GeV/c^. The three charged tracks are required to form a vertex. 

We select B — > Kx'y —>■ K*7i'y candidates (here and throughout this section, K* denotes 
-ft'* (892) for simplicity) by requiring the invariant mass of K^7c~ to be within ±75 MeV/c^ 
of the nominal K* mass. The resulting Mbc and AE distributions are shown in Figs. ^ and 
H, respectively. Using the same fitting procedure as is used for the B -^ i^2 (1430)7 analysis, 
we obtain 46.4 ± 7.3 lll'^ events from the Mbc distribution. This is consistent with the yield 
obtained from the AE distribution. These events are dominated by 5 — *> K*ti'-) as a K* 
mass peak is clearly seen in the Kir invariant mass distribution in Fig. |^ However, other 
contributions, which can arise either from B^ — > K^ p^'^ or non-resonant B^ — > K^'k^t^^'-j, 
are not negligible. We estimate these contributions to be 5.8 ±2.2 IqJ events from the region 
1.1 GeV/c^ < Mk-k < 1.4 GeV/c^. An additional contribution from other b ^ s'y decay is 
estimated to be 0.9 ±0.6 events from MC. After subtracting these non K* contributions, we 
obtain a B+ ^ K*^Tr+-f yield of 39.7 ± 7.4 1^;^ events. 

From the Kx invariant mass [M^x) distribution (Fig. ^), we observe a broad structure 
below 2.0 GeV/c^ that can be explained, for example, as a sum of two known resonances 
around 1.4 GeV/c^ and 1.7 GeV/c^, but cannot be explained by a single known resonance 
or phase space decay. We observe no excess above 2.0 GeV/c^, indicating that the Mkx < 
2.0 GeV/c^ cut does not introduce a significant inefficiency. 

To estimate the efficiency oi B ^ K*ti'-), we analyze B -^ Ki( 1400)7 ^ud B — >■ 
ii'*( 1680)7 MC samples, use the mean of the efficiencies as the central value, and assign 
the difference to the systematic error. As a result, the efficiency becomes (3.13 ± 0.47)% 
including the other systematic errors in Table |I[ We determine the B -^ K*7T'~f branching 
fraction, 

B{B -^ K*7rr, Mk*^ < 2.0 GeY/c^) = (5.6 ± 1.1 ± 0.9) x 10"^ 

There are four known resonances, 7^1(1270), 7^1(1400), K*{U10) and K^{U30), that can 
contribute to the signal around Mkx = 1-4 GeV/c^. In the region of 1.2 GeV/c^ < M^^ < 
1.6 GeV/c^, we obtain 22.9 ± 5.1 +!;? events from the Mbc distribution. The ir2*(1430) 
contribution is estimated to be 2.6 ± 1.4 events from our branching fraction measurement, 
assuming B^ — > i^2(1430)°7 and B^ -^ 7^2(1430)^7 have equal branching fractions. The 
reconstruction efficiencies are about the same for Ki(1400)7 and i^*(1410)7 and a factor 
of two lower for Ki(1270)7 including the iri(1270) -^ Kp and Ki(1270) -^ ir*(1430)7r 
contributions. We interpret the signal yield as an upper limit on the weighted sum of the 
three resonances, 

^B{B -^ Ki{1270)-f) + B{B -^ is:i(1400)7) ± B{B -^ K*(1410)7) 

< 5.1 X 10"^ (90% C.L.) 

This limit on B{B -^ i^*(1410)7) is more stringent than that obtained from the cos6'hei 
distribution of Krcy decays. 
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Next, we select B -^ Kx'y — > Kp'y candidates by requiring the invariant mass of the 7r"''7r~ 
combination to be within ±250 MeV/c^ of the nominal p mass. To veto B -^ Kx'y —^ K*7T'y 
events, we reject a candidate if the invariant K^tt~ mass is within ±125 MeV/c^ of the 
nominal K* mass. The Mbc distribution and the Kx invariant mass distribution are shown 
in Figs. P and |l^, respectively. From the Mbc distribution, we obtain a signal yield of 
24.5 ± 6.4 ^2!! events. We subtract the contribution of 2.3 ± 1.2 events from other 6 — > 57 
decays. 

The Mkx spectrum of these events (Fig. |TD|) shows a large peak around 1.7 GeV/c^. 
One possible explanation is a large i^*( 1680)7 component with a small contribution from 
i^i(1270), but since there are quite a few resonances around 1.7 GeV/c^, a detailed analysis 
will be required to disentangle the resonant substructure. The reconstruction efficiency for 
B -^ Kp'j, which is M^^ dependent, is determined to be (1.51 ± 0.25)% by assuming a 
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candidates. 
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mixture of 7^1(1270) and K*{1Q80) with a ratio from the Mkx fit result. So far we find 
no signal outside the p mass window; neglecting the non-resonant Kinx'-y contribution, we 
determine the B — > Kp'-y branching fraction, 

B{B ^ Kpr, Mkp < 2.0 GeV/c^) = (6.5 ± 1.7 1\\) x 10"^ 

The i^p7 final state in the mass range around 1.3 GeV/c^ is effective for the search 
of i? — >• i^i(1270)7, because the i^i(1270) has a large ((42 ± 6)%) branching fraction to 
Kp. The second largest contribution is from i^2(1430) -^ Kp ((8.7 ± 0.8)%) from which 
we expect less than one event. Other contributions are much smaller. We search for B — »• 
ii'i(1270)7 decays by requiring \Mkx — ^i^i(i27o)| < 0.1 GeV/c^, as shown in Fig. 0. We 
find 4 candidates in the signal box with a background expectation of 1.19 events. Using a 
reconstruction efficiency of (0.41±0.06)%, we obtain an upper limit oiB{B — > ^^1(1270)7) < 
9.6 X 10-5 (90% C.L.). 



V. CONCLUSION 



We have searched for radiative B meson decays into kaonic resonances that decay into a 
two-body or three-body final states together with a high energy photon. We observe sizable 
signals in 5 — ;> ii'g (1430)7, B — > K*7r'j and B — > Kp'j decays and determine the branching 
fractions for these channels. The measured branching fractions respectively correspond to 
about 4%, 17% and 19% of the total 6 -^ 57 branching fraction assuming the SM calculation 
IP or existing measurements p,|10|,|ll[ as the denominator. Adding 15% from the K*{892)^ 
branching fractions, these decay modes sum up to about half of the entire b —>■ s'j process. 

For the Kn'-f final state, the i^2 (1^30)7 component is separated from a possible 
i^* (1410)7 or non-resonant contribution using a helicity angle analysis. 

For the three-body final states, we observe B —>■ K*7['~f and B -^ Rpj signals separately 
for the ffist time; however, the possible contribution of many kaonic resonances prevents 
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us from further identification of such resonances with the current statistics. We find no 
significant signal for B -^ 7^1(1270)7 decay in the Kp'-f final state. 
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